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We have revisited the fragment-asperity interaction model recently introduced by Sotolongo-Costa 
and Posadas (Physical Review Letters 92, 048501, 2004) Q by considering a different definition for 
mean values in the context of Tsallis nonextensive statistics and introducing a new scale between the 
earthquake energy and the size of fragment ear 3 . The energy distribution function (EDF) deduced 
in our approach is considerably different from the one obtained in the above reference. We have 
also tested the viability of this new EDF with data from two different catalogs (in three different 
areas), namely, NEIC and Bulletin Seismic of the Revista Brasileira de Geofi'sica. Although both 
approaches provide very similar values for the nonextensive parameter q, other physical quantities, 
e.g., the energy density differs considerably, by several orders of magnitude. 
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I. INTRODUCTION 

Over the last two decades, a great deal of atten- 
tion has been paid to the so-called nonextensive Tsallis 
entropy, both from theoretical and observational view- 
points. This particular nonextensive formulation 0, Q 
seems to present a consistent theoretical tool to investi- 
gate complex systems in their nonequilibrium stationary 
states, systems with multifractal and self-similar struc- 
tures, systems dominated by long-range interactions, 
anamolous diffusion phenomena, among others. Some 
recent applications of Tsallis entropy S q ^\ to a number 
of complex scenarios is now providing a more definite 
picture of the kind of physical problems to which this 
g-formalism can in fact be applied. 

In this regard, systems of interest in geophysics has 
also been studied in light of this nonextensive formal- 
ism. In this particular context, the very first investiga- 
tion was done by Abe Q who showed that the statistical 
properties of three-dimensional distance between succes- 
sive earthquakes follow a g-exponential function with the 
nonextensive parameter lying in the interval [0,1] 
Since then, other geophysical analyses have been per- 
formed as, for instance, the statistics of the calm time, 
which indicates a scale-free nature for earthquake phe- 
nomena and corresponds to a ^-exponential distribution 
with q > 1 [(|, and models for temperature distribu- 
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tions and radon emission of volcanos |jj • More recently, a 
very interesting model for earthquakes dynamics related 
to Tsallis nonentensive framework has been proposed by 
Sotolongo-Costa and Posadas (SCP Model) Q. Such a 
model consists basically of two rough profiles interact- 
ing via fragments filling the gap between them, where 
the fragments are produced by local breakage of the lo- 
cal plates. By using the nonextensive formalism the au- 
thors of Ref. not only showed the influence of the size 
distribution of fragments on the energy distribution of 
earthquakes but also deduced a new energy distribution 
function (EDF), which gives the well-known Gutenberg- 
Richter law Q as a particular case. 

However, in dealing with this nonextensive framework, 
a particular attention must be paid to the possible def- 
initions for mean values, which play a fundamental role 
within the domain of this nonextensive statistics . In 
this concern, recent studies of the properties of the rela- 
tive entropy and the Shore- Johnson theorem for consis- 
tent minimum cross-entropy principle, revealed the ne- 
cessity of the so-called g-expectation value in studies in- 
volving this nonextensive statistical mechanics (see [Tol ] 
for details). Thus, by introducing this g-definition of 
mean value we re-analyzed the fragment-asperity inter- 
action model of Sotolongo-Costa and Posadas [lj. More- 
over, a new scale law between the released relative energy 
e and the 3-dimensional size of fragments has also been 
introduced. By using the standard method of entropy 
maximization we also deduced a new energy distribu- 
tion function, which differs considerably from the one 
obtained in Ref. Q. In order to test the viability of 
our appoach we used data taken from two seismic cata- 
logs, namely, NEIC and Bulletin Seismic of the Revista 
Brasileira de Geofisica. It is shown that although both 
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approaches provide very similar values for the nonexten- 
sive parameter q, the other physical quantity, e.g., the 
energy density differ by several orders of magnitude. 

This paper is organized as follows. In Sec. II, the stan- 
dard formalism of nonextensive statistical mechanics is 
reexamined, as well as the theorical basis of SCP model. 
In Sec III, a new EDF is analytically calculated through 
extremization of Tsallis' entropy under the constrains of 
the ^-expectation value and normalization condition. In 
Sec IV, we test this new EDF with data from two dif- 
ferent catalogs and estimate the best-fit values for the 
nonextensive parameter q and the proportionality con- 
stant between the released relative energy e and the vol- 
ume of the fragments r 3 , i.e., the energy density, a. We 
end this paper by emphasizing the main results in the 
conclusion Section. 



II. NON-EXTENSIVE FRAMEWORK AND SCP 
MODEL 



also necessary since the process of violent fractioning is 
very probably a nonextensive phenomenon, leading to 
long-range interactions among the parts of the object 
being fragmented (see, e.g., PI El). In reality, such an 
influence was earlier emphasized in other investigations 
[Tij . In general lines, the SCP model follows similar ar- 
guments to those presented in Refs. 0] being, however, 
a more realistic seismic model than the one proposed in 
Ref. 0| . In particular, the theoretical ingredients reads: 

• the mechanism of relative displacement of fault 
plates is the main cause of earthquakes; 

• the surfaces of the tectonic plates are irregular and 
the fragments filling the space between them are 
very diverse and have irregular shapes; 

• the mechanism of triggering earthquakes is estab- 
lished through the combination of irregularities of 
the fault planes and the distribution of fragments 
between them; 



In this Section, we recall the nonextensive theoretical 
basis of the SCP model. As widely known, the Tsal- 
lis' statistics generalizes the Botzmann-Gibbs statistics 
in what concerns the concept of entropy. Such formalism 
is based on the parametric class of entropies given by 



S q= £i = -k B J ' P q {<y) \n q p(<j)da, 



(1) 



where ks is the Boltzmann constant. In the SCP model, 
p(<y) stands for the probability of finding a fragment of 
relative surface a (which is defined as a characteristic 
surface of the system), q is the nonextensive parameter 
and the q-logarithmic function above is defined by 



\n q p = (1 - g )-l(pi-* - 1), (P>0) 



(2) 



which recovers the standard Boltzmann-Gibbs entropy 
Si = — fc_B f plrtpd 3 p in the limit q — > 1. It is worth 
mentioning that most of the experimental evidence sup- 
porting Tsallis proposal are related to the power-law dis- 
tribution associated with S q ^\ descripition of the classi- 
cal TV-body problem 

The SCP model is a simple approach for earthquakes 
dynamics revealing a very interesting application of the 
Tsallis' framework. Indeed, the fundamental idea con- 
sists in the fact that the space between faults is filled 
with the residues of the breakage of the tectonic plates. 
In this regard, the authors studied the influence of the 
size distribution of fragments on the energy distribution 
of earthquakes. The theoretical motivation follows from 
the fragmentation phenomena [T^ | in the context of the 
geophysics systems. In this latter work, Englaman et al 
showed that the standard Botzmann-Gibbs formalism, 
although useful, cannot account for an important feature 
of fragmentation process, i.e., the presence of scaling in 
the size distribution of fragments, which is one of the 
main ingredients of the SCP approach. Thus, a nonex- 
tensive formalism is not only justified in SCP model but 



• the fragment distribution function and conse- 
quently the EDF emerges naturally from a nonex- 
tensive framework. 

From the above arguments, the EDF deduced in Ref. 
is given by 



1"S(-V ,„) = loj- A' ( ^— | I 



(3) 



x log 1 + a(q - 1)(2 - q)VTv x 10 2 



According to Ref. the above expression describes 
very well the energy distribution in all detectable range of 
magnitudes, unlike the empirical formula of Gutenberg- 
Richter @. 



III. NEW APPROACH 

Now, let us discuss the standard method of maximiza- 
tion of the Tsallis entropy. Here and hereafter, the Boltz- 
mann constant is set equal to unity for the sake of sim- 
plicity. Thus, the functional entropy to be maximized is 

SS* q =s(s q + aj p{(j)da - pa^J = 0, (4) 

where a and (3 are the Lagrange multipliers. The con- 
strains used above are the normalization of the distribu- 
tion 



p{a)da = 1 



and the q-expectation value 



<J q =< a > q = I aP q (a)d(j 



(5) 



(6) 
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FIG. 1: The relative cumulative number of earthquakes [Eq. 1121 1 as a function of the magnitude m. In all panels, the data 
points correspond to earthquakes lying in the interval 3 < m < 8. (a): Samambaia fault - Brazil: 100 events from Bulletin 
Seismic of the Revista Brasileira de Geoh'sica. (b): New Madrid fault - USA: 173 data points taken from NEIC catalog, (c): 
Anatolian fault - Turkey: 8980 events from NEIC catalog. The best-fit values for the parameters q and a are shown in the 
respective Panel. A summary of this analysis is also shown in Table I. 



with the escort distribution |17| given by 

P ^ (7) 

By considering the same physical arguments of Ref. , 
we derive, after some algebra, the following expresion for 
the fragment size distribution function 



(1-g) 
(2-9) 



(a -<j q 



(8) 



which corresponds to the area distribution for the frag- 
ments of the fault plates. Here, however, differently from 
Ref. [lj, which assumes e ~ r, we use a new energy scale 
e ~ r 3 . Thus, the proportionality between the released 
relative energy e and r 3 (r is the size of fragments) is now 
given by a — a q = (e/a) 2 ^ 3 , where a scales with r 2 and a 
(the proportionality constant between e and r 3 ) has di- 
mension of volumetric energy density. In particular, this 
new scale is in accordance with the standard theory of 
seismic rupture, the well-known seismic moment scaling 
with rupture length (see, for instance 0). 

The new EDF of earthquakes is, therefore, obtained by 
changing variables from a — a q to (e/a) 2 ^ 3 . From JSJ it 



is straightforward to show that 

Ce~ids 



p(e)de = 



(9) 



[1 + C'e 2 /3] 

which has also a power-law form with C and C" given by 



C = 



3a 2 / 3 



and C = -- 



(1-9) 



(10) 



(2- g )a 2 / 3 ' 

In the above expression, the energy probability is written 
as p(e) = 71(e) /N, where n(e) corresponds to the num- 
ber of earthquakes with energy e and TV total number 
earthquakes. 



IV. TESTING THE NEW EDF WITH THE 
CUMULATIVE NUMBER OF EARTHQUAKES 

In order to test the viability of the new EDF above 
derived [Eq. JJJJl] we introduce the cumulative number of 
earthquakes, given by integral 1] 



N 



p(e)de, 



(11) 
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TABLE I 


: Limits to q 


and a 




Fault 


Ref. 


q 


a 


California - USA 




1.65 


5.73 x 10~ 6 


Iberian Penisula - Spain 


1.64 


3.37 x 10~ 6 


Andalucfa - Spain 


m 


1.60 


3.0 x 10 -6 


Samambaia - Brazil 


This Paper 


1.60 


1.3 x 10 10 


New Madrid - USA 


This Paper 


1.63 


1.2 x 10 10 


Anatolian - Turkey 


This Paper 


1.71 


2.8 x 10 10 



where N E> is the number of earthquakes with energy 
larger than e. Now, substituting into Ijllfl . and con- 
sidering m — j log e (m stands for magnitude) it is pos- 
sible to calculate the above expression. In reality, note 
that depending on the value of q the limits of the integral 
(|ll|l presents a cutoff on the maximum value allowed for 
energy e, which is given by e max = v / a 2 / 3 (2 - q)/(l - q) 
for the intervals q < 1 and q > 2, while for 1 < q < 2 
the cutoff is absent in the distribution. Note also that in 
the limit q — ► 1, e rnax — > oo and p(e) goes to the expo- 
nential function. As matter of fact, the calculation of the 
integral (|11|) for g ^ 1 leads to the general expression 



log(AT >m ) = \ogN 



xlog 



2-g 

i-q 

1- g 

2- g 



(12) 



10 2 



,2/3 



which, similarly to the modified Gutenberg-Ricther law 
(See, e.g., Refs. for more details), describes appropri- 
ately the energy distribution in a wider detectable range 
of magnitudes. 

Figure 1 shows the relative cumulative number of 
earthquakes (G m > = N m> /N) as a function of the mag- 
nitude m. The data points, corresponding to earthquakes 
events lying in the interval 3 < m < 8, were taken from 
two different catalogs, namely, Bulletin Seismic of the 
Revista Brasileira de Geofi'sica (left panels) and NEIC 
(central and right panels). The left, central and right 
Panels show the results of our analysis for the Samam- 
baia fault, Brazil (100 events), New madrid fault, USA 
(173 events), and Anatolian fault, Turkey (8980 events), 
respectively. We note that, similarly to original version 
of SCP model, our approach, represented by Eqs. 



(|12fl . provide a very good fit to the experimental data 
of the two catalogs here considered. It is worth empha- 
sizing, however, that the energy density differ by several 
orders of magnitude from our model to the original SCP 
model. Therefore, we expect that other independent es- 
timates of the parameter a may indicate which approach 
is more physically realistic. The estimates of the param- 
eters q and a obtained in this paper and in Ref. are 
summarized in Table 1. 

V. CONCLUSION 

In Ref. ^3 ' what seems to be the correct definition for 
expectation values within the Tsallis noncxtcnsivc statis- 
tical mechanics was rediscussed. Based on properties of 
the generalized relative entropies and the Shore-Johnson 
theorem, it was shown that the expectation value of any 
physical quantity in this extended framework converges 
to the normalized ^-expectation value, instead of to the 
ordinary definition. 

In this paper, by considering this necessity of q- 
expoctation values in Tsallis nonextensive framework, we 
have revisited the fragment-asperity interaction model 
for earthquakes, as introduced in Ref. 0. A new energy 
distribution function has been calculated, which allowed 
us to determine the relative cumulative number of earth- 
quakes as a function of the magnitude. Additionally, a 
new scale law between the released relative energy e and 
the volume of fragments r 3 has also been introduced, i.e., 
in agreement with the so-called seismic moment scaling 
with rupture length. As discussed earlier^ although our 
analysis and the one presented in Ref. [lj provide very 
similar values for the nonextensive parameter q, the other 
physical quantity, e.g., the energy density differ by several 
orders of magnitude. It would be interesting, therefore, 
if we could have experimental estimates for these quan- 
tities in order to compare the predictions of the models. 
Finally, it is worth mentioning that the estimates for the 
nonextensive parameter from the two catalogs here con- 
sidered (Fig. 1) are consistent with the upper limit q < 2, 
obtained from several independent studies involving the 
Tsallis nonextensive framework |20|. 
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